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We experimentally study the fluid flow induced by a broad, penetrable barrier moving through 
an elongated dilute gaseous Bose-Einstein condensate. The barrier is created by a laser beam swept 
through the condensate, and the resulting dipole potential can be either attractive or repulsive. We 
examine both cases and find regimes of stable and unstable fluid flow: At slow speeds of the barrier, 
the fluid flow is stationary due to the superfluidity of the condensate. At intermediate speeds, we 
observe a non- stationary regime in which the condensate gets filled with dark solitons. At faster 
speeds, soliton formation completely ceases and a remarkable absence of excitation in the condensate 
is seen again. 
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The fluid flow past an obstacle is one of the most proto- 
typical experiments studying superfluidity. Accordingly, 
moving an obstacle through a superfluid has been met 
with considerable interest, both in the context of su- 
perfluid helium [1] as well as in the context of dilute 
gaseous Bose-Einstein condensates (BECs) [2]. In ex- 
periments with BECs, the role of the obstacle can be 
played by a laser beam that creates a dipole potential 
for the atoms. By moving a small, strongly repulsive 
dipole beam through a BEC, evidence for a critical veloc- 
ity above which superfluidity breaks down was obtained 
^, 4]. In those experiments, the laser was impenetrable 
for the atoms, and the diameter of the laser beam was 
chosen smaller than the size of the BEC. The atoms could 
therefore flow past the sides of the moving beam. The 
critical velocity observed in the experiments was much 
lower than that given by the well-known Landau criterion 
[5]. Vortices can be produced by such narrow obstacles, 
and indeed they were observed experimentally [6] , and in 
numerical simulations [7]. Likewise, in the context of op- 
tical lattices, heating due to dissipative motion has been 
used as a tool to study superfluidity in the presence of 
periodic potentials [8]. 

In this Letter we consider a situation that is comple- 
mentary to the experiments described in 0, 0]: In our 
case, the BEC is contained in a narrow, elongated trap 
and a barrier much wider than the radial extent of the 
BEC is moved through it (Fig. [TJa)). The atoms cannot 
flow past the sides of the barrier, and we use a penetrable 
barrier allowing atoms to flow through the barrier. The 
barrier is created by the dipole potential of an ellipti- 
cally shaped laser beam that is swept along the long axis 
of the elongated BEC. Such a geometry approximates a 
one-dimensional problem, and it is theoretically expected 
that in this case solitons play a similar role to the role 
played by vortices in higher dimensions [§, 11, EH- 
Recently, shedding of solitons has also been discussed in 
terms of Cerenkov-radiation [13|, [l^, ll5|. In this Let- 
ter we provide experimental evidence for the existence of 



several different flow regimes in the wake of the moving 
barrier. For slow and for fast velocities, the wake behind 
the moving barrier appears unperturbed. In an interme- 
diate velocity regime solitons appear. This behavior is 
observed with repulsive barriers as well as with attrac- 
tive potentials. 

The starting point for the experiments described in 
this paper are elongated ^''Rb BECs of about 4.5 • 
10^ atoms in the \F = I^tuf = —1) state, held in 
a loffe-Pritchard type trap with trapping frequencies 
{u;x/{'^'7r)^u;yz/{27r)} = {7.1,203}Hz. The axis of weak 
confinement (x-axis, see Fig. [U^a)) is oriented horizon- 
tally. Evaporative cooling is performed until no thermal 
cloud of uncondensed atoms is visible. A laser that is far 
detuned from the Rb absorption lines is directed along 
the second horizontal axis (y-axis). Depending on the 
detuning, it can create an attractive or repulsive dipole 
potential for the atoms. This dipole beam can be moved 
along the x-axis at various speeds, inducing flow in the 
BEC. The dipole beam waist in the vertical z-direction 
is much larger than the radial extent of the BEC so that 
the intensity variation of the beam along the z-direction 
can be neglected. The imaging direction coincides with 
the direction of the dipole beam. A notch filter prevents 
the dipole beam from being detected by the camera. 

As we show in this Letter, solitons can be produced 
by sweeping the dipole beam through the condensate. 
The width of these solitons is below our spatial imaging 
resolution when the BEC is held in the trap. Therefore 
we employ a fast antitrapping technique to expand the 
BEC before imaging, similar to the procedure described 
in [16]. The expansion imaging procedure starts with a 
microwave adiabatic rapid passage, converting the atoms 
from the \F = l^mp = —1) state to the antitrapped \F — 
2^mF = —2) state. We then set the magnetic bias field 
to 48 G and let the cloud expand for 3 ms. During this 
expansion, the magnetic field is positioned such that the 
antitrapping potential keeps the cloud from falling under 
the influence of gravity. The atoms are then imaged in a 
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FIG. 1: (a) Axis convention. The dipole beam propagates in 
the -y direction and is moved through the BEG along the -x 
direction, (b) Onset of excitation at slow sweep speeds, calcu- 
lated from images with parameters as in (c)-(k). A repulsive 
dipole beam is swept through the cloud from the right to the 
left. The white arrow in Fig. (c) indicates the end position. 
Sweep speeds are (c) 0.4 mm/s, (d) 0.5 mm/s, (e) 0.6 mm/s, 
(f) 0.7 mm/s, (g) 0.8 mm/s, (h) 1 mm/s, (i) 1.3 mm/s, 
(j) 2 mm/s, (k) 3.3 mm/s. Images show a 200 jj.m wide region. 
Other parameters see text. 



100 G bias field using resonant absorption imaging. The 
antitrapping predominantly leads to a radial expansion, 
so that the aspect ratio of the BEG changes from Ryz/Rx 
= 29 for the trapped BEG to a value of about 1.6 at the 
end of the expansion time. 

In a first set of experiments, shown in Fig. [H we study 
the flow induced by a weak repulsive potential swept 
through the BEG at various speeds. The potential was 
produced by a laser beam with a wavelength of 660 nm 
and a power of 37 jiiW. The beam shape was elliptical 
with waists Wx = 7.6 /im and Wz = 25.8 jam. Such a 
laser produces a weak repulsive potential whose height is 
about 24% of the chemical potential of the condensate. 
The healing length in the center of our typical conden- 
sates is 0.17 jam so that the dipole beam can be consid- 
ered wide with respect to the healing length even in the 
direction of the sweep. The dipole beam starts out on 
the right, outside the BEG, and is then swept towards 
the left at a constant speed as given in the figure cap- 
tion. The sweep is stopped at the position indicated by 
the white arrow in Fig.[Hc), and the BEG is imaged after 
antitrapped expansion as described above. 

We can identify several different regimes, depending on 



the sweep speed. For very slow sweep speeds, such as in 
Fig. [TJc), the part of the BEG through which the dipole 
beam has been swept appears essentially unaffected, and 
the only visible effect of the beam is the density suppres- 
sion that it leaves at its end position. Since the BEG 
is superfluid, it is not surprising that it remains unaf- 
fected by a slowly moving penetrable barrier. At faster 
sweep speeds (Fig. [Hd)-(g)), dark notches appear in the 
wake of the dipole beam. We interpret these notches as 
the solitons predicted in theoretical studies of the one- 
dimensional problem d, [13, El, E, [HI (see below for 
an experimental verification). Some of the notches ob- 
served in the wake of the dipole beam show slight bend- 
ing or nonuniform contrast along the plane, while others 
are very straight and uniform. In the darkest solitons, 
we observe a suppression of the central density by about 
50%. However, it is conceivable that the true depth of 
some of the solitons is larger than this because slight 
misalignments of the imaging direction with respect to 
the soliton planes can reduce the apparent depth. In- 
terestingly, fewer solitons are observed when the sweep 
speed is increased beyond 0.9 mm/s (Fig. [TJh). This ef- 
fect is accompanied by the appearance of an increasingly 
wider region of low density directly behind the dipole 
beam (Fig. [Hh),(i)). At faster sweep speeds, the con- 
densate appears again unaffected by the obstacle, as seen 
in Fig. [IJk). The apparent absence of excitations in 
the wake of the beam is very striking. Suppression of 
phonons and solitons at hig h obstacle velocities was the- 
oretically discussed in jiol. [ill. In particular, [ll| 
describes an explanation of related numerical findings 
by using an analogy to the radiation of capillary-gravity 
waves in a classical fluid. 

To quantify the amount of excitation in the BEG, we 
first calculate the smooth, overall shape of each BEG 
by considering only the lowest Fourier components of 
each image. Then the root-mean-square deviation of each 
original image from the smooth, overall shape is used as a 
measure of the excitation present in the BEG. To reduce 
the influence of the inhomogeneous density distribution 
along the x-axis, we only consider a region in the center of 
the BEG that is 35 /am wide in the x-direction and 70 /am 
wide along the z-direction. The resulting quantity, ^, is 
plotted versus the barrier speed in Fig.fT^b). The points 
at zero velocity are taken from unperturbed BEGs in the 
absence of a dipole beam, and mostly indicate imaging 
noise. The onset of excitations above a critical velocity 
of roughly 0.3 mm/s can clearly be seen in this plot. We 
do not extend this plot to velocities above 0.8 mm/s be- 
cause at higher speeds the broad low-density region that 
develops in the wake of the barrier (as seen in Fig. [Tfh), 
(i)) affects our measure. For comparison, the bulk speed 
of sound in the center of our BEGs is 3 mm/s. In a 
cigar shaped BEG, averaging over the tightly confined 
direction leads to a lower speed of sound [13, [3, 0. 
In our case, this results in a speed of 2.1 mm/s which 
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FIG. 2: Images (upper row) and corresponding integrated 
cross sections (lower row) for various endpoints of the sweep 
as indicated by the white arrows. Laserpower 69 /xVK, sweep 
velocity 1.4 mm/s. 



is still significantly larger than the observed critical ve- 
locity. This situation is reminiscent of the observation 
of a surprisingly low critical velocity in vortex shedding 
experiments in [3|, U]. For a more precise analysis, the 
complicated fluid flow through the barrier would have to 
be taken into account. In addition, the nonuniform ra- 
dial density profile of the BEG is expected to lead to a 
deviation from the Bogoliubov dispersion relation and to 
a lower critical velocity, as discussed in [21|. 

When a more strongly repulsive potential is used, we 
observe similar behaviors. However, solitons appear al- 
ready at lower sweep velocities, and the velocity above 
which the flow becomes stationary again increases. Also, 
a more pronounced accumulation of atoms in front of 
the dipole beam is seen, delimited by a steep wavefront 
(Fig. [2]). For weaker potentials the number of solitons 
in the BEG is reduced. For example, when using dipole 
beams with an intensity of 25 /iVK, often only one or two 
individual solitons are produced. 

In a second set of experiments we now show that the 
dark notches in the wake of the dipole beam are not just 
soundwaves but do indeed behave like solitons. To prove 
this, we compare the stability of the observed notches 
with the stability of a simple density suppression in the 
BEG (Fig. [3]). To obtain a very clearly visible den- 
sity suppression, a relatively strong dipole beam with a 
power of 129 /nW was swept through half of the BEG. 
The height of the dipole potential was about 85% of the 
chemical potential of the condensate so that a deep den- 
sity suppression was left in the BEG at the end position 
of the sweep. For clarity this suppression is marked by 
the white arrow in Fig. [31(a). The sweep velocity was 
adjusted to 0.2 mm/s, leading to pronounced soliton for- 
mation in the wake of the laser. At the end of the sweep, 
the dipole beam was switched off, and the BEG was al- 
lowed to evolve in trap for a variable evolution time be- 
fore starting the expansion imaging procedure. Typical 
images for various evolution times are shown in Fig. [3l 
At very short evolution times (Fig. [31(a)), the right part 
of the cloud, through which the dipole beam has been 
swept, shows stacks of dark notches, and the density sup- 




FIG. 3: Stability of the observed solitons. Times for which the 
BEG is allowed to evolve in trap after turning the dipole beam 
off at the end of the sweep and before starting the expansion 
imaging procedure: (a) ms, (b) 10 ms, (c) 20 ms, (d) 30 ms, 
(e) 40 ms, (f) 50 ms. The arrow in (a) indicates the end 
position of the sweep. Other parameters see text. Images 
show a 200 /xm wide region. 



pression produced by the dipole beam at its end position 
is clearly visible. This density suppression gets filled in 
by the BEG over a time scale of about 10 ms (Fig. [31(b)). 
Gonsidering that the longitudinal speed of sound in our 
cigar-shaped BEG is 2.1 mm/s, it is very plausible that 
the density suppression with an initial width of roughly 
20 fim closes during 10 ms. However, the much narrower 
dark notches in the right part of the cloud remain visible 
for much longer times and slowly spread out (Fig. [3Kb)- 
(f)). After evolution times of over 50 ms they are still 
visible and have spread over the entire BEG (Fig. [S^f)). 
We often observe very straight dark notches at this time. 
Then their contrast gradually fades and they disappear. 
These measurements demonstrate the stability of the pro- 
duced dark notches, suggesting that they are indeed very 
distinct from sound waves and correspond to the solitons 
predicted by theory for the one-dimensional problem. 

Our condensates are elongated, but not one- 
dimensional. Therefore, the question arises which influ- 
ence the transverse dimension has on the existence of 
solitons. In previous experiments, individual solitons in 
a BEG with repulsive interactions have been engineered 
using phase imprinting |22|, |23| or wavefunction engineer- 
ing [24] . They have also been observed in the evolution of 
a density gap in a BEG produced by a stopped light tech- 
nique [25]. One important outcome of the phase imprint- 
ing studies is that the exact nature of the solitary texture 
depends on the strength of the radial trapping. The nu- 
merical studies given in [2^ reveal that if the chemical 
potential exceeds ten times the level spacing in the tightly 
confined direction, i.e. if jn/hcOyz > 10, phase imprinting 
produces a nonstationary kink that decays into phonons 
on a time scale of 1/ujyz. If this ratio is below five, the 
kink decays into a proper dark soliton. Experimentally, 
solitons produced by phase imprinting were found in ex- 
periments with ji/hwyz = 7 [22]. While our mechanism 
for soliton formation is not based on phase imprinting, 
a similar influence of the transverse confinement on the 
existence of solitons can be expected. In particular, for 
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our BECs fi/hujyz = 9.7, which is not too different from 
the experiments in [22[. 

By changing the laser wavelength to 830 nm we can 
also study the effects of an attractive potential, in which 
case the density is increased instead of decreased in the 
region of the potential. As in the case of repulsive poten- 
tials, theories predict soliton formation at intermediate 
speeds and vanishing excitation at very slow and very 
fast speeds jlo|- We indeed observe such regimes. How- 
ever, deviating from our observations with repulsive po- 
tentials, we have not found a region of suppressed density 
behind an attractive dipole laser that would correspond 
to Fig. [TJi). Results obtained with an attractive dipole 
beam are shown in Fig. [H A dipole beam with a wave- 
length of 830 nm and a waist of 8.2 /im in the x-direction 
and 31.7 jam in the z-direction was used. For the analysis 
in Fig. [41(a), a laser power of 10 jiiW was chosen, leading 
to a potential depth of about 17% of the chemical po- 
tential of the BEC. As before we quantify the excitation 
by plotting ^ versus the sweep speeds. Since we do not 
observe a low-density region corresponding to Fig. HJi), 
we can now use ^ as a measure for the excitation over 
the whole range of velocities studied. Despite the weak- 
ness of the perturbing potential, soliton formation again 
occurs at surprisingly low velocities. 

For a laser power of 10 fiW, only a few solitons are 
produced, as shown in Fig. [41(b). If the laser power is in- 
creased to 19 /iW^, the soliton formation can become quite 
pronounced (Fig. [41(c). This demonstrates that not only 
a repulsive potential, but also an attractive potential can 
be an efficient tool to produce solitons. One-dimensional 
theories for the case of attractive potentials predict the 
onset of soliton formation to occur right at the speed of 
sound ^o[. However, we see that already at a speed of 
1.25 mm/s, below the central speed of sound, the central 
part of the BEC can be filled with solitons. This dis- 
crepancy may hint at an influence of the inhomogeneous 
density across the condensate in our experiments. 

In conclusion, we have observed stationary and non- 
stationary fluid flow regimes when a penetrable barrier 
is moved through a BEC. A slowly moving barrier leads 
to stationary flow. At intermediate velocities, the flow is 
non-stationary and is marked by the emergence of soli- 
tons in the wake of the barrier. At fast velocities, the 
soliton shedding stops and the flow appears stationary 
again. In a one-dimensional geometry, the solitons play 
a role analogous to vortices present in flow around obsta- 
cles in higher dimensions [3,, 4, . Thus our experi- 
ments complement these previous studies that have been 
hallmark experiments for investigating the superfluidity 
in dilute gaseous BECs. 
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FIG. 4: Effect of an attractive potential swept through the 
condensate, (a) Excitation strength versus sweep speed for a 
laser power of 10 /xVK, other parameters see text, (b) Laser 
power 10 /iiW, sweep speed 1.25 mm/s. (c) Laser power 
19 /iVK, sweep speed 1.25 mm/s. Images show a 200 /xm wide 
region. The white arrow in (b) indicates the end position of 
the sweep. 
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